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Remote-sensing  reflectance  t Rn),  which  is  defined  as  the  ratio  of  water-leaving  radiance  (Lw)  to  down- 
welling  irradiance  just  above  the  surface  (^(0+)),  varies  with  both  water  constituents  (including  bottom 
properties  of  optically -shallow  waters)  and  angular  geometry.  Lw  is  commonly  measured  in  the  field  or  by 
satellite  sensors  at  convenient  angles,  while  Erf(0+)  can  be  measured  in  the  field  or  estimated  based  on 
atmospheric  properties.  To  isolate  the  variations  of  Rn  (or  Lw )  resulting  from  a  change  of  water  consti¬ 
tuents,  the  angular  effects  of  Rn  (or  Lw)  need  to  be  removed.  This  is  also  a  necessity  for  the  calibration 
and  validation  of  satellite  ocean  color  measurements.  To  reach  this  objective,  for  optically-deep  waters 
where  bottom  contribution  is  negligible,  we  present  a  system  centered  on  water’s  inherent  optical  proper¬ 
ties  (IOPs).  It  can  be  used  to  derive  IOPs  from  angular  R n  and  offers  an  alternative  to  the  system 
centered  on  the  concentration  of  chlorophyll.  This  system  is  applicable  to  oceanic  and  coastal  waters 
as  well  as  to  multiband  and  hyperspectral  sensors.  This  IOP-centered  system  is  applied  to  both  numeri¬ 
cally  simulated  data  and  in  situ  measurements  to  test  and  evaluate  its  performance  The  good  results 
obtained  suggest  that  the  system  can  be  applied  to  angular  Rn  to  retrieve  IOPs  and  to  remove  the 
angular  variation  of  Rn.  ©  2011  Optical  Society  of  America 
OCIS  codes:  010.4450,  010.0280. 


1 .  Introduction 

Water-leaving  radiance  (LW)  Wm_2sr_1.  See  Table  1 
for  primary  symbols  used  in  this  article)  is  the  light 
intensity  per  solid  angle  resulting  from  absorption 
and  scattering  processing  in  the  water.  Spectral  Lw 
provides  key  information  to  remotely  sense  water 
constituents  in  the  upper  water  column,  as  it  varies, 
both  in  magnitude  and  spectral  shape,  with  changing 
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concentrations  of  optically  significant  constituents 
and  bottom  properties  (bottom  depth  and  bottom 
reflectivity)  if  the  sea  bottom  is  optically  shallow. 
Accurate  determination  of  Lwy  from  in  situ  or  satel¬ 
lite  sensors,  is  a  fundamental  goal  for  ocean  color 
satellite  missions. 

Lw  varies  not  only  with  wavelength  U,  nm)  and 
water  constituents  but  also  with  angular  geometry: 
0Si  <P  (collectively  represented  as  O  in  the  follow¬ 
ing  text  unless  further  clarification  is  needed).  Here, 
0S  is  the  solar  zenith  angle  in  air  ( 6V ,  <p),  determines 
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Table  1.  Notation 


Symbol  Definition  or  Description 

Units 

a 

Absorption  coefficient 

m"1 

b 

Scattering  coefficient 

m  1 

b„ 

Backscattering  coefficient  (=  bbw  +  bbp) 

m  1 

b but 
bbp 

Backscattering  coefficient  of  water  molecules 
and  particles,  respectively 

m  1 

bf 

Forward  scattering  coefficient 

m‘ 1 

c 

Beam  attenuation  coefficient  (=  a  +  b) 

m  1 

Dd 

Distribution  function  of  downwelling  irradiance 

Ed 

Downwelling  irradiance 

Wm'2 

Eod 

Downwelling  scalar  irradiance 

Wnr2 

Eou 

Upwelling  scalar  irradiance 

Wm-2 

f 

Model  coefficient  for  subsurface  irradiance 
reflectance 

- 

a 

Model  coefficients  for  subsurface  remote 
sensing  reflectance 

s  r 1 

G 

Model  coefficients  for  above-surface  remote 
sensing  reflectance 

sr-1 

Attenuation  coefficient  of  upwelling  radiance 

m  1 

L 

Radiance 

W  m"2  sr-1 

K 

Upwelling  radiance 

Wm'2  sr"1 

Water-leaving  radiance 

W  m~2  sr"1 

[I'U'W 

Normalized  water-leaving  radiance 

W  m-2  sr-1 

n 

Refractive  index  of  water 

- 

Q 

Ratio  of  irradiance  to  radiance 

sr 

rn 

Subsurface  remote  sensing  reflectance 
(=Lu(0-)/£d(0-)) 

sr' 1 

Rr. 

Above-surface  remote-sensing  reflectance 

(—  Lw(0+)/Ed(0+)) 

sr'1 

P 

Volume  scattering  function 

m'1  sr'1 

K 

Quasi-diffuse  attenuation  coefficient  (—  a  +  bb) 

m  1 

X 

Wavelength 

nm 

Os 

Solar  zenith  angle  in  air 

rad 

0v 

Sensor  nadir-view  angle  in  air 

rad 

<P 

Sensor  azimuth  angle  in  relation  to 
the  solar  plane 

rad 

n 

Collectively  representing  sun-sensor 
angular  geometry  (0S  ,6V  q>) 

rad 

p 

Fresnel  reflectance 

- 

!K 

Divergence  factor  when  radiance  enters  into  air 
from  below  the  surface 

the  angular  direction  of  Lw  with  0V  the  zenith  angle 
of  the  radiance  (=  sensor's  view  angle  from  nadir)  in 
air,  and  q>  the  azimuth  angle  relative  to  the  sun 
plane.  Figure  1  presents  the  angular  coordinates 
used  in  this  article.  For  ft  defined  by  0S  =  30°,  0V  = 
30°  and  q>  =  180°,  Lw  is  in  a  direction  going  to  the  sun 
(the  sensor  has  its  back  face  towards  the  sun  and 
measures  more  backscattered  photons  from  the 
sun  beam);  while  for  an  ft  with  0S  =  30°,  0V  =  30°, 
and  q>  —  0°,  Lw  is  in  a  direction  going  further  away 
from  the  sun  (and  a  remote  sensor  would  measure 
the  specular  reflection  of  the  solar  light  from  air- 
sea  interface). 

Because  LW(X ,  ft)  varies  with  ft  for  the  same  water 
constituents  under  the  same  solar  illumination  [1],  it 
is  necessary  to  remove/correct  this  angular  variation 
(also  called  the  bidirectional  reflectance  distribution 
function  [BRDF])  in  order  to  isolate  the  constituent- 
caused,  versus  geometric  (different  time  or  different 
location)  variation  of LW(X ,  ft).  This  is  also  important/ 


Fig.  1.  (Color  online)  Angular  coordinate  used  in  this  study  and 
the  IOP-centered  system. 


required  for  the  vicarious  calibration  of  space  sensors 
and  the  validation  of  Lw(X,£l)  derived  from  different 
sensors  which  have  different  measurement  geome¬ 
tries.  Because  of  this  demand,  the  concept  of  normal¬ 
ized  water-leaving  radiance  [LW(X)]N  was  introduced 
[2-4].  This  quantity  is  a  hypothetical  value  of  Lw 
with  zenith  angle  of  0°  (sensor  looking  vertically 
down,  0V  =  0°),  the  sun  at  zenith  (0s  =  0°),  and  a 
black  sky  (no  atmosphere  between  the  sun  and  the 
water).  Changes  in  [LW(X)]N  at  different  locations 
and  times  thus  represent  only  changes  of  water  con¬ 
stituents.  Because  the  actually  measured  property  is 
LW(X,Q.)>  a  method  to  convert  Lw(Xy£l)  to  [L^i)]^  is 
required. 

In  the  1990s,  for  optically  deep  waters,  Morel  et  al. 
[1,4,5]  developed  a  conversion  system  for  optically 
deep  waters  based  on  the  Case  1  assumption  [6], 
which  uses  concentration  of  chlorophyll  (Chi)  as  an 
input  for  calculations.  To  briefly  summarize,  a 
relationship  between  [L^i)]^  and  LW(X,Q)  is  ex¬ 
pressed  as 

[Z/wMW  _  FqW _ 

LW(X,  ft)  Ed( 0M,  0S)  m(0s ,  0V)  [Q(X)}Nf(X ,  os)  ‘ 

(i) 

Here,  F0  (W  m~2)  is  the  solar  irradiance  at  the  top  of 
the  atmosphere  and  for  the  mean  sun-earth  distance, 
and  Ed  is  the  spectral  downwelling  irradiance  just 
above  the  surface  when  Lw  is  measured.  Parameter 
91  represents  cross-surface  effects  [5,7]  when  upwel- 
ling  radiance  enters  into  air  from  the  water;  f  relates 
subsurface  irradiance  reflectance  with  the  ratio  of 
backscattering  coefficient  (6^,  m"1)  to  absorption 
coefficient  (a,  m_1)  [8];  and  Q  represents  the  ratio 
of  subsurface  upwelling  irradiance  to  subsurface 
upwelling  radiance  [1]  corresponding  to  LW(X,CL ). 
Consistently,  [fU)]^  and  [Q(i)]jv  indicate  the  values 
of  f(X, 0S)  and  Q(i,ft)  corresponding  to  [L^A)]^. 

In  Eq.  (1),  F0  is  known,  Ed  can  be  calculated  based 
on  known  atmosphere  information  [9],  variation  of  91 
can  be  determined  based  on  information  of  angular 
geometry  and  is  found  nearly  independent  of  wind 
speed  [7].  The  remaining  unknowns  are  f  and  Q  as 
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they  vary  with  wavelength,  angular  geometry,  and 
water  constituents.  Based  on  the  Case  1  assumption 
[10,11],  where  the  water's  spectral  optical  properties 
can  be  determined  by  the  concentration  of  Chi  alone 
[  12] ,  variations  of  both  f  and  Q  (and  their  ratios)  have 
been  developed  [5]  and  validated  [13]  for  a  set  of  Chi 
values  and  wavelengths.  Consequently,  when  Chi  is 
known  or  estimated  from  measured 
[L^yl)]^  can  then  be  determined  based  on  Eq.  (1). 

This  system  requires  knowledge  of  the  concentra¬ 
tion  of  Chi  and  works  only  with  Case  1  waters 
[1,5,14],  thus  its  application  excludes  waters  that 
do  not  fall  into  the  Case  1  category,  which  includes 
many  coastal  as  well  as  oceanic  waters  [15].  The 
following  details  an  alternative  approach  that  uses 
the  water's  inherent  optical  properties  (IOPs)  as  in¬ 
puts  for  calculations  and  does  not  require  the  Case  1 
assumption,  thus  it  could  be  applied  to  both  oceanic 
waters  and  coastal  waters.  The  performance  of  this 
approach  to  both  numerically  simulated  data  and 
field  measured  data  is  also  presented.  Because  bot¬ 
tom  reflectance  may  significantly  modify  the  BRDF 
of LW(X,  O)  [16],  the  study  here  is  focused  on  optically- 
deep  waters.  In  addition,  the  contribution  to  Lw(X,£l) 
resulted  from  inelastic  scattering  (like  Raman  scat¬ 
tering  [17,18]  or  chlorophyll-a  fluorescence  whose 
angular  distribution  is  isotropic  [19])  is  not  included. 

2.  Theoretical  Basis 

The  subsurface  remote-sensing  reflectance  (r^,  sr-1) 
is  defined  as 


rn{W) 


Lui  o-,A,n') 

Ed(0-,X,ffsY 


(2) 


where  Lu(0“,A,fi')  is  the  subsurface  upwelling  radi¬ 
ance  and  Ed( O~,A,0^)  is  the  subsurface  downwelling 
irradiance  ( ffs  is  tfie  in-water  zenith  angle  of  9S). 
LW(X,Q)  can  then  be  expressed  as  [20] 

Lw{X,tt)  =  1  ~  y)  gd(0-,  x,  ffs)rn(X,  no.  (3) 

where  p  is  the  Fresnel  reflectance  corresponding  to  ffv 
(in-water  zenith  angle  of  0V)>  and  n  is  the  refractive 
index  of  water.  This  relationship  indicates  that  the 
angular  variation  of  LW(X,Q)  is  mainly  determined 
by  the  angular  shape  of  rre(i,n'). 

Based  on  the  radiative  transfer  theory,  Zaneveld 
[21,22]  has  shown  that  r  „,(>!,  O')  of  optically-deep 
waters  is 


_ _ _ 

-c(X)+kL(W)-fL(W)bf(X) 

fo*  fo/2P(n'>  d)L(X,  O')  sin(0')d0'd<p' 

£orf(0-^//s) 


Here,  Dd  is  the  distribution  function  for  the  downwel¬ 
ling  light  [23],  c  (m-1)  is  the  beam  attenuation  coeffi¬ 
cient,  ki  (m-1)  is  the  diffuse  attenuation  coefficient 


for  upwelling  radiance  [22],  fL  is  the  shape  param¬ 
eter  for  radiance  [22],  bf  (m_1)  is  the  forward 
scattering  coefficient,  ft  (nr*  sr_1)  is  the  volume  scat¬ 
tering  function  (VSF)  of  the  water  media,  L  is  radi¬ 
ance  at  an  angle  (0',$'),  and  E ^  is  the  downwelling 
scalar  irradiance.  In  essence,  rre  depends  on  the 
water  IOPs  and  the  radiance-sun  angular  geometry 
[1,22,24]. 

Under  the  quasi-single  scattering  assumption  [25], 
Eq.  (4)  for  zenith-going  radiance  (Lu(0~,  A,  0°))  can  be 
simplified  as  [26,27] 


rre(A,  0°) 


Dd(X,ff8)fi(180-ffs) 
a(X)+kL(X,  0°) 


(5) 


This  formulation,  although  with  larger  errors 
resulting  from  the  omission  of  multiple  scattering, 
highlights  that  rre(>l,  O')  is  fundamentally  dependent 
on  the  optical  properties  [5,24]  (e.g.,  absorption  and 
backscattering  coefficients)  and  the  angular  shape 
of  the  VSF  (or  phase  function)  [27,28].  Consequently, 
it  is  necessary  to  know  both  in  order  to  accurately 
remove  the  angular  variation  of  rre(A,  O')  or  LW(X,  H), 
especially  for  wavelengths  where  single  scattering 
dominates. 

Unfortunately,  in  ocean  color  remote  sensing,  the 
only  available  information  is  Lu  (A,H)  (when  bound¬ 
ary  conditions  are  adequately  known),  therefore  the 
optical  properties  of  the  water  must  be  derived  first 
from  in  order  to  achieve  the  BRDF  correc¬ 

tion.  While  the  angular  shape  of  p  is  composed  of  both 
water  molecule  and  particle  scattering,  it  has  been 
found  that  the  angular  shape  of  the  particle  scatter¬ 
ing  phase  function  is  not  a  constant  and  varies  with 
particle  size  and  composition  [29—31],  To  extend  the 
dilemma  further:  at  present,  there  is  no  reliable 
method  yet  to  accurately  estimate  the  particle  phase 
function  shape  from  remote-sensing  measurements. 
As  a  result,  in  ocean  color  remote  sensing,  a  phase 
function  (or  phase  functions  [5])  for  particle  scatter¬ 
ing  has  to  be  assumed  in  order  to  develop  an  analy¬ 
tical  model  [1,23,24,32]  to  process  remotely  sensed 
data.  Fortunately,  the  variation  of  particle  phase 
function  of  global  waters  is  within  limited  ranges 
[31].  The  fact  that  particle  phase  functions  do  vary 
indicates  that  the  BRDF  correction  in  ocean  color  re¬ 
mote  sensing  will  always  be  a  best-effort  estimate, 
unless  multiple  scattering  has  created  a  completely 
diffuse  upwelling  light  field. 

As  discussed  above,  even  after  we  assume  that  the 
particle  scattering  phase  function  is  known,  it  is  still 
required  to  know  the  optical  properties  (e.g.,  absorp¬ 
tion  and  diffuse  attenuation  coefficients)  for  the 
correction  of  angular  effects.  Because  the  diffuse  at¬ 
tenuation  coefficients  are  mainly  determined  by  the 
absorption  (a)  and  backscattering  ( bb )  coefficients 
[33-35],  thus  the  logical  focus  is  to  derive  a  and  b ^ 
from  the  measured  LW{X,  O).  This  could  not  be 
achieved  through  Eq.  (4)  as  it  is  not  directly  inverti¬ 
ble.  To  expedite  the  inversion  process,  numerous 
attempts  have  been  made  to  obtain  a  simplified, 
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but  sufficiently  accurate,  analytical  formulation  to 
describe 

3.  Model  for  Remote-Sensing  Reflectance 

In  the  past  decades,  following  various  numerical  si¬ 
mulations,  a  series  of  formulas  have  been  developed 
to  achieve  a  simplified  expression  for  Eq.  (4).  These 
include  the  pioneer  polynomial  model  developed  by 
Gordon  et  aL  [24]  for  nadir- viewed  r 


rrsM=0°)  = 


go  +£i 


(  Mj)  Y  _ 


bbW 


(6) 


Albert  and  Mobley  [36]  expanded  the  above  to  a 
4th-order  polynomial,  with  coefficients  developed 
for  various  angles  and  wind  speed  ( w ,  m  s_1): 


<7) 


To  explicitly  account  for  the  effects  of  molecular 
and  particle  scatterings  [37],  Lee  et  al.  [38,39]  pro¬ 
posed  a  two-term  formula  to  model 


rrsOLH')  =gw(nf) 


bbwW 

a(A)  +  bb{X) 


+  gP(A,LY) 


t>bPW 

a(A)  +  bb(X) 


(8) 


Separately,  Park  and  Ruddick  [40]  used  a  4th-order 
polynomial  with  model  coefficients  as  functions  of  the 
relative  contributions  (vb)  of  molecular  scattering  to 
the  total  scattering: 


^(-1,(1')  =  Jftffl'.Ifc) 


i= 1 


bb(A) 


a(A)  +  bb(X) 


(9) 


Lately,  Van  Der  Woerd  and  Pasterkamp  [41] 
proposed  a  significantly  different  polynomial  form 


For  an  easier  description  and  following  the  ap¬ 
proach  of  Gordon  et  aL  [24],  R „  of  optically  deep 
waters  is  considered  as  a  generic  function  of  a  model 
parameter  (G)  and  water’s  IOPs: 

R„(i l,  fi)  =  G(i l,  fi)Fun(IOP(A)).  (12) 

Here  the  angular  variation  of  R is  represented 
by  G,  which  may  also  vary  with  IOPs  as  shown  for 
subsurface  remote-sensing  reflectance  (see  Eq.  (6)). 
The  remote-sensing  reflectance  corresponding  to 

[Lw}n  is 


[RMn  =  [GWW  Fun(IOPW),  (13) 

and  [L^^)]^  will  simply  be  the  product  of  [RrsWW 
and  F0(A).  The  challenge  of  the  BRDF  correction  then 
becomes  how  to  convert  G(A,  fl)  to  [G(A)]N  when 
/{n(A,  fl)  is  the  only  known,  where  as  G(A,  D)  may  also 
vary  with  IOPs  (see  Eqs.  (6H9)). 

With  an  objective  to  efficiently  process  a  large  vol¬ 
ume  of  satellite  data,  and  at  the  same  time  explicitly 
reveal  the  different  effects  of  molecular  scattering 
versus  particle  scattering  [37,39],  we  selected  a  for¬ 
mula  similar  as  Eq.  (8)  to  model  H): 

fl)  =  (g^A)  +  Gf  (SI) 

+  (GS(n,  +  G?(n,^H)^,  (14) 

with  K  =  a  +  bb.  This  formulation  integrates  the 
cross-surface  effects  (31(0$,  0J)  and  the  expression 
for  /Vs  into  one  equation,  simply  to  expedite  data 
processing.  This  is  also  because  the  value  of  31(0$,  0V) 
can  be  well  determined  based  on  angular  geometry 
and  it  is  nearly  independent  of  the  wind  speed  [7]. 
Expressions  for  R^/LO)  with  explicitly  separated 
$R(0s,0v)  and  rre,  however,  are  also  available. 


ln[rre(A,n'))  =  ££*></(  <nM«W)]'MW))F- 

1  J  1 

(10) 

These  models,  although  not  exhaustive  of  those  in 
the  literature,  all  used  the  particle  phase  function 
measured  by  Petzold  [29]  and  in  general  produce 
consistent  r for  given  IOPs.  They  do,  however,  have 
subtle  differences  in  complexity  and  representation 
of  radiative  transfer. 

Because  a  remote  sensor  measures  L^(>l,n)  in¬ 
stead  of  Lu(0“,/i,n'),  we  take  a  slightly  different 
approach  in  the  modeling  step.  Defining  the  re¬ 
mote-sensing  reflectance  (R^)  as  the  ratio  of  LW(X ,  D) 
to£d(OV,05): 


Lw  (a,  n) 

Ed(0\X,0s)' 


(11) 


Fig.  2.  Particle  phased  function  shapes  (Petzold  average  and  1% 
Fournier  and  Forand)  used  for  the  numerical  simulation  The 
open  diamond  represents  one  weighted  average  of  the  particle 
assemblage. 
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Table  2.  Sample  Values  of  (G^fl),  Gf  ffl),  Gj(n),  Gf  (H))  for  Angular  Rn 


0s  => 

0° 

0V  =*> 

0° 

<p=> 

0° 

G° 

0  0604 

Gl 

0.0406 

G°p 

0,0402 

G\ 

0.1310 

0° 

15° 

0 

O 

CO 

o 

O 

CO 

90° 

90° 

0.0596 

0.0590 

0.0516 

0.0562 

0.0408 

0,0411 

0.1420 

0.1461 

o 

O 

CO 

0° 

o 

o 

CO 

o 

O 

Tt 

90° 

135° 

0.0584 

0,0581 

0.0601 

0.0581 

0.0418 

0  0414 

0.1492 

0.1458 

15° 

p 

o 

CO 

© 

O 

o 

O 

135° 

135° 

00614 

00624 

0.0524 

0.0524 

0.0425 

0.0434 

0.1408 

0,1406 

The  model  coefficients  (Gq  (H),  Gq(O),  and 

G?(fl))  in  Eq.  (14)  are  now  dependent  on  angular 
geometry  and  phase  function,  but  independent  of 
absorption  and  backscattering  coefficients  or  wave¬ 
length.  Similarly  to  earlier  studies  [8,23,32,42],  the 
values  of  (Gj^D),  G^'(fl),  G^fD),  and  Gj(fi))  were 
derived  with  numerically  simulated  remote-sensing 
reflectance  and  with  preselected  particle  phase 
functions. 

We  took  advantage  of  the  synthetic  dataset  in¬ 
cluded  in  the  IOCCG  Report  #5  [43],  where  500  sets 
of  spectral  IOPs  are  available  (bb/(a  +65):  0.0002- 
0.37)  in  a  spectral  range  from  400  to  800  nm  (resolu¬ 
tion  is  10  nm).  Two  particle  phase  functions  were 
used  for  the  simulation,  one  for  mineral  like  particles 
and  the  other  for  phytoplankton,  to  better  represent 
the  collective  particle  scattering  of  an  assemblage  of 
particles.  The  former  is  represented  with  the  Petzold 
average  [20],  while  the  latter  is  represented  by  a 
particle  phase  function  with  a  backscattering- 
to-scattering  ratio  as  1%  [44].  Figure  2  presents 
the  Petzold  average,  the  1%  Foumier-Forand  phase 
function,  and  an  example  of  the  blended  phase  func¬ 
tions.  The  simulation  (by  Hydrolight  [45])  in  the 
IOCCG  Report  #5  [43]  covered  only  two  sun  angles 
(30°  and  60°);  here,  we  extended  that  to  cover  sun 
angles  from  0  to  75°. 

For  0S  from  0  to  75°  (resolution  is  15°),  0V  from  0  to 
70°  (resolution  is  10°),  and  (p  from  0  to  180°  (resolu- 


Fig.  3.  Schematic  data  flow  chart  of  the  IOP-centered  BRDF 
correction  scheme. 


tion  is  15°),  a  total  (6x7x13-4  6)  set  of  constants 
representing  (G^  (O),  G^(n),  Gg(n),  and  Gj(H))  were 
derived  (by  least-square  fit)  from  the  simulated 
i?re(A,  fl).  Table  2  provides  a  few  samples  of  the  model 
coefficients.  Note  that  these  values  can  be  easily  re¬ 
vised/updated  when  a  more  suitable  phase  function 
(s)  is  accepted  by  the  community. 

4.  Retrieve  of  iOPs  from  Angular  Rn  and  the  BRDF 
Correction 

As  discussed  in  Section  3  and  based  on  Eq.  (14),  the 
absorption  and  backscattering  coefficients  are  re¬ 
quired  in  order  to  obtain  [7?,^]#.  Figure  3  shows  a 
general  data  flow  chart  of  the  IOP-based  BRDF  cor¬ 
rection  scheme.  We  here  take  the  quasi-analytical  al¬ 
gorithm  (QAA)  strategy  [46,47]  for  the  derivation  of 
the  IOPs  from  R^A,  H).  This  is  because  in  the  QAA 
system,  the  R^A,  O)  modeled  with  the  derived  a  and 
bb  represents  exactly  the  measured  R^A,!!).  On 
the  other  hand,  if  a  spectral  optimization  approach 
[48-50J  is  used  for  the  derivation  of  a  and  bby  the 
modeled  R^A,!!)  with  the  derived  a  and  bb  does 
not  match  the  measured  (A,  D)  exactly.  This  differ¬ 
ence  in  Rre(A,n)  will  result  in  extra  uncertainty 
in  the  estimated  [RrsWW  when  the  optimization- 
derived  a  and  bb  is  applied. 

In  the  QAA,  absorption  at  a  reference  wavelength 
Ao,  a(Ao),  is  estimated  first,  and  then  particle  back- 
scattering  at  this  wavelength  (&6n(Ao)).  a  (A)  is  then 
derived  analytically  from  R^A,  O)  after  propagating 
this  6&p0*o))  to  bbp(X)  with  a  power-law  model  [10,46]. 
These  derived  a  (A)  and  bbp(X)  are  introduced  into 
Eq.  (14),  and  applying  the  (GjJ'(fl),  G^(H),  Gg(fi), 
and  G^(H))  values  for  the  angular  geometry 
(0°,0°,0°)  of  normalized  water-leaving  radiance  to 
determine  [R^A)]^. 

In  the  derivation  process,  A0  is  selected  as  550,  555, 
or  560  nm  for  Moderate  Resolution  Imaging  Spectro- 
radiometer  (MODIS),  Sea-viewing  Wide  Field-of- 
view  Sensor  (SeaWiFS),  and  MEdium  Resolution 
Imaging  Spectrometer  (MERIS)  sensors,  respec¬ 
tively,  and  a(Ao)  is  estimated  from  (QAA-v5,  http:// 
www.  ioccg.  org/gro  ups/Soft ware_OC  A/Q  ;AA_v5 .  pd  0 


a(/lo)  =  aw(lQ)  +  10-1  146-1366,-0  469/ 


(15) 


The  coefficients  in  Eq.  (15)  are  averages  for  those 
sensors,  but  individual  coefficients  for  each  sensor 
are  available.  Because  the  variation  of  the  spectral 
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(21) 


ratio  of  R^A,  D)  for  different  D  is  quite  limited,  the 
constants  in  Eq.  (15)  are  used  for  R n  of  all  D. 
Impact  of  the  error  of  a(Ao)  on  the  BRDF  correction 
is  presented  in  Subsection  5.B.3. 

When  a(Ao)  is  known,  Eq.  (14)  becomes  a  quadratic 
function  of  bbp(A0): 

Ax2+Bx  +  C  =  0 ,  (16) 

with  x  =  bbp(A 0)  and 


(A  =  GP  +  G?-Rn, 

{  B  =  Gq  bbw  4-  G%(a  4-  bbw)  -  2Rrs(a  +  bbw), 

{  C  =  G%bbw(a  +  bhw)  -Rn(a  +  bbwf  + 

(17) 

bbp(Ao)  is  then  derived  by  solving  Eq.  (16).  Note 
that  the  Rre,  a,  and  bbw  in  Eq.  (17)  are  values  at 
A$.  bbp(A)  is  further  derived  with  a  power-law  function 
[10]  where  the  exponent  is  estimated  empirically  [46] 
(also  see  http://www.ioccg.org/groups/Software_OCA/ 


QAA_v5.pdf). 

When  bbp(A)  is  known,  Eq.  (14)  becomes  a 
tic  equation  of  k: 

quadra- 

RrsK2  -Xk-Y  =  0, 

(18) 

with 

jX  =  G-bhw  +  Gp0bb 
\Y  =  G‘?(bbw)*  +  G{(bbpf. 

(19) 

And 

Q  ^  bbw  b bp , 

(20) 

after  k  is  derived  by  solving  Eq.  (18).  bbw ,  bbp ,  and 

k  in  Eqs.  (18M20)  are  values  at  A. 


8  = 


\l-lw(a,o°)\ 
LW{A,  0°) 


x  100. 


Here,  L  represents  either  LW(A,  Q)  (before  correction) 
or  (A,  Degenerated  Lw(A,0°)  (i.e.,  after  correction). 
The  distribution  of  6  is  shown  in  Fig.  4. 

It  is  found  that  without  a  BRDF  correction,  only 
23.7%  of  the  dataset  have  <5  within  5%,  and  more  than 
50%  of  the  data  have  8  >  10%,  which  is  consistent 
with  results  shown  in  earlier  studies  [1,5].  After 
applying  the  IOP-based  BRDF  correction  scheme, 
however,  92.3%  of  the  data  have  <5  within  5%,  and 
only  1.4%  of  the  data  have  8  >  10%.  These  results 
demonstrate  that  the  IOP-based  BRDF  correction 
scheme  (which  is  composed  of  the  Rn  model  along 
with  the  IOP  retrieval  algorithm)  works  very  well 
with  the  numerically  simulated  data.  This  could  be 
in  part  because  the  model  and  dataset  have  the  same 
particle  phase  functions. 

B  Field  Measured  Data 


1.  Comparison  of  NuRADS  Measurements 

The  upwelling  radiance  distribution  was  measured 
using  two  systems.  For  the  AOPEX  dataset  (experi¬ 
ments  in  the  Mediterranean  Sea,  August  2004), 
NuRADS  [51]  were  used.  For  the  Monterey  Bay  ex¬ 
periment,  a  combination  of  the  NuRADS  systems 
were  used,  which  also  allowed  measurement  of  the 
polarized  radiance  distribution  [52].  The  NuRADS 
systems  measure  the  in-water  upwelling  spectral 
radiance  distribution  and  are  based  on  a  fisheye 
camera  system  and  an  electronic  camera.  With  these 
systems,  the  entire  hemispherical  radiance  distribu¬ 
tion  can  be  measured,  at  one  wavelength,  in  a  single 
image  [51]. 

For  the  data  used  in  this  paper,  the  upwelling 
radiance  distribution  was  averaged  over  10  min 


5.  Validation  of  the  lOP-Based  BRDF  Correction 
Scheme 

The  IOP-based  BRDF  correction  scheme  is  applied  to 
both  numerically  simulated  and  field  measured  data 
to  evaluate  its  performance. 

A.  Data  Simulated  with  Hydrolight 

For  a  self-consistency  check,  the  above  BRDF  correc¬ 
tion  scheme  was  first  applied  to  a  subset  of  data  si¬ 
mulated  by  Hydrolight  (Section  3).  We  focused  on  Lw 
with  the  sun  at  60°  from  zenith,  and  the  sensor  view¬ 
ing  angle  in  a  range  of  30-70°  and  azimuth  angles  of 
0-180°,  an  angular  domain  which  shows  a  large  dif¬ 
ference  between  LW(A,  0°)  (nadir-measured  Lw)  and 
LW(A,  fl).  Here,  we  compared  L^^H)  with  LW(A,  0°) 
(i.e.,  if  there  is  no  correction  at  all),  and  compared 
them  between  (A,  H) -generated  LW(A,  0°)  with 
LW(A ,  0°)  (i.e.,  after  a  BRDF  correction).  For  each  pair 
of  LWJ  an  absolute  difference  is  calculated  as 


- 
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Fig.  4.  Distribution  of  absolute  percentage  difference  (<5)  of  Lw. 
Simulated  data  were  used,  and  focused  on  the  angular  domain 
that  6S  —  60°,  6V  —  30°-70°,  and  <p  —  0o-~180o.  Solid  circle, 
no  BRDF  correction;  open  diamond,  after  IOP-based  BRDF 
correction. 
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periods.  Each  NuRADS  image  takes  on  the  order  of 
0.1s  to  acquire;  however,  because  of  the  time  re¬ 
quired  to  rotate  through  the  spectral  bands  and  store 
the  data,  there  are  approximately  five  images  at 
each  wavelength  which  can  be  averaged  in  the 
10  min  period. 

While  LU(X,  D)  is  measured  in  radiometric  units, 
because  of  the  measurement  technique  (all  angles 
in  a  single  image),  the  angular  shape  of  the  upwelling 
radiance  (LU(X,  Cl)/Lu(X,  0°))  is  a  much  more  accurate 
parameter.  In  this  paper  we  compare  the  angular 
shape  of  Lw(X,Cl)/Lw(X,  0°)  in  the  model  and  mea¬ 
surements.  For  measurements,  LW(X,  n)/Lw(X,0°) 
was  obtained  by  applying  the  ratio  of  5R(fl)//?(0°) 
(see  Gordon  [7])  to  Lu(X,Cl)/Lu(X,  0°),  which  was  cal¬ 
culated  directly  from  the  NuRADS  measurements.  In 
the  model,  LW(X,Q)/LW{X,  0°)  was  estimated  from 
/2rs(>l,  Q)  using  the  IOP-centered  approach. 

Results  of  two  completely  different  waters  are 
presented,  where  the  sky  was  clear  without  clouds. 
One  was  measured  in  blue  water  (Mediterranean 
Sea,  AOPEX  station  0807125900,  Aug.  7,  2004)  with 
the  sun  at  ~30°  from  zenith;  the  other  was  measured 
in  green  water  (Monterey  Bay,  Sept.  15,  2006)  with  a 
sun  angle  ~60°  from  zenith.  For  the  AOPEX  station, 
Rn  was  derived  from  upwelling  radiance  and  down- 
welling  irradiance  measured  from  three  TriOS  radio¬ 
meters  (for  example,  see  Lubac  and  Loisel  [53]).  For 
the  Monterey  Bay  station,  was  measured  at 
(30°,  90°)  with  a  handheld  spectroradiometer  [54,55]. 
The  measured  of  the  two  stations  are  presented 
in  Fig.  5.  The  absorption  coefficient  at  440  nm  (de¬ 
rived  from  QAA)  of  the  two  stations  was  ~0.025  m-1 
and  0.9  m-1,  while  their  euphotic  zone  depths  [56] 
were  ~108m  and  7.2  m,  respectively,  indicating 
significantly  different  water  environments. 

Following  the  IOP-centered  BRDF  correction 
scheme,  the  absorption  and  backscattering  coeffi¬ 
cients  were  derived  first  from  by  incorporating 
their  corresponding  model  coefficients  for  (Gq^D), 
G?  (n),  G£(H),  and  G?(n)).  With  these  IOPs,  R^X,  0°) 
andi?rs(>l,  D),  and  the  ratio  fy/R^X,  0°)  (which 
equals  to  Lw(X,n) /Lw(X,0°))  were  calculated  with  the 


angular  coefficients  developed  in  Section  3.  As  exam¬ 
ples,  Figs.  6  and  7  (blue  symbols)  compare  measured 
and  calculated  LW(X,  0.)/Lw(X ,  0°)  of  selected  viewing 
angles  and  wavelengths. 

For  the  AOPEX  station  (Fig.  6),  because  the  sun  is 
at  30°  from  zenith,  the  angular  variation  among 
LW(X,  n)/Lw(X,0°)  is  quite  narrow  (between  0.9  and 
1.1)  for  0V  within  70°.  In  this  relatively  clear  water, 
molecule  scattering  dominates  the  scattering  pro¬ 
cesses  in  the  shorter  wavelengths  (b^w/b^  ~  75%  at 
412  nm),  and  the  distribution  of  upwelling  radiance 
is  more  isotropic  [1].  In  addition,  the  smaller  solar 
zenith  angle  limits  the  contribution  of  forward  scat¬ 
tered  photons  from  the  solar  beam  to  L^(A,60°,0°). 

To  evaluate  the  overall  performance  of  the 
BRDF  correction  scheme  in  the  remote  sensing 
domain  (0V  =  0-71°,  (p  —  0-180°),  as  examples, 
Lw(X,Os,Gv,(p)/Lw(X,0s,Oo,Oo)  for  X  as  412  and 
550  nm  are  presented  at  the  lower  panel  of  Fig.  6,  re¬ 
spectively.  Statistically,  the  coefficient  of  determina¬ 
tion  (R2)  is  0.86  and  0.77,  with  averaged  S  values  of 
1.4%  and  1.7%  for  the  two  wavelengths,  respectively. 

For  comparison,  Lw(Xyti)/Lw(X, 0°)  ratios  derived 
from  the  Chl-based  approach  [4,5]  are  also  presented 
(green  square  in  Fig.  6),  where  Chi  concentration  was 
~0.08mg/m3  based  on  the  OC4v4  algorithm  [57]. 
The  same  R(il)/R( 0°)  ratios  used  for  the  NuRADS 
data  were  applied  to  convert  the  below-surface 
radiance  ratio  to  the  above-surface  radiance  ratio. 
Generally,  the  results  of  the  two  approaches  are  simi¬ 
lar,  which  is  evidenced  by  the  R2  values  (0.89  and 
0.84  for  412  nm  and  550  nm,  respectively)  and  the 
S  values  (1.4%  and  2.2%  for  412 nm  and  550 nm, 
respectively).  These  results  imply  that  for  oceanic 
waters,  a  consistent  BRDF  correction  could  be 
achieved  from  both  the  Chl-based  and  IOPs-based 
approaches,  but  more  tests  and  evaluations  with 
extensive  measurements  are  necessary. 

For  the  measurements  made  in  Monterey  Bay, 
however,  a  different  picture  emerges.  For  the 
measurements  with  0^  =  30°,  the  1^(440, 30°,  <p)/ 
L^(440,0°,0°)  ratio  is  generally  around  1.0,  except 
for  those  where  cp  is  close  to  0°,  an  angle  where  more 
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Fig.  5.  Remote  sensing  reflectance  of  the  two  sample  stations:  (left)  measured  in  the  Mediterranean  Sea,  (right)  measured  in  the 
Monterey  Bay. 
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forward  scattered  solar  photons  contribute  to  Lw. 
This  is  significantly  pronounced  for  0V  =  60°  and  (p 
close  to  0°  (Fig.  7).  Because  the  sun  was  at  -60° 
(-40°  below  the  sea  surface),  and  because  particle 
scattering  dominates  in  this  environment  (b^/b^ 
10%  at  412  nm),  more  multiple  forward  scattered  so¬ 
lar  photons  contribute  to  LW(X,  60°,  q?  -  0°)  than  to 
L„(i,0°,0°),  and  consequently  LW(X,  60°,  g?  -  0°)/ 
Lwy,0°,0°)  was  as  high  as  1.7  [58].  Note  that,  for 
0V  =  60°  and  <p-90°,  Lw(;t,60o,p - 0°)/L„(A,0°,0°) 
is  only  around  1.2.  For  0V  =  60°  and  cp  -  180°,  how¬ 
ever,  because  it  is  in  the  angular  range  of  the 
instrument  self  shadow,  the  measured  LW(X,  60°, 
60°,  <p-  180°)/!^  (A,  60°,  0°,0°)  dropped  to  -0.80. 
Instrument  self  shading  [59],  however,  is  minor  in 
the  blue  spectral  region  for  clear  waters  [51]. 

Similarly  to  Fig.  6,  LW(X,  0s,  0VJq>)/Lw(A, 0s,  0°,  0°) 
for  X  as  440  and  550  nm  are  presented  at  the  lower 
panel  of  Fig.  7,  respectively  The  R 2  values  are 
0.95  and  0.98,  with  averaged  S  values  as  3.5%  and 
4.1%  for  the  two  wavelengths,  respectively.  Because 
of  instrument  self  shading,  data  with  (p  in  the  range 
of  150-180°  (the  open  symbol  in  Fig.  7)  were  excluded 
in  the  statistical  analysis.  The  Chi  value  was 
-27.7  mg/m3  (derived  empirically  from  the  ra¬ 
tio),  which  is  beyond  the  Chi  range  of  the  Chl-based 
approach  [5],  thus  no  LW(X,  0S,  0v,q>)/Lw(X,  Gs,  0°,  0°) 
ratio  was  derived  from  the  Chl-based  approach  for 
this  measurement. 


The  above  results  indicate  that  between  the  NuR- 
ADS  measurements  and  the  R ^  estimations,  the 
Lw(6s,Ov,<p)/Lw(0s,Oo,O°)  ratios  of  both  oceanic 
and  coastal  waters  were  very  consistent  qualitatively 
and  quantitatively  and  thus  validate  the  IOP- 
based  BRDF  correction  system,  at  least  for  these 
measurements. 

2.  Comparison  Using  the  Wire-Stabilized 
Environmental  Profiling  Radiometer  and  the  Ocean 
Color  Component  of  the  Aerosol  Robotic  Network 
Data 

Another  evaluation  was  carried  out  with  63  pairs  of 
water-leaving  radiance  data  determined  from  almost 
simultaneous  above-  and  in-water  radiometric  mea¬ 
surements  (i.e.,  collected  within  15  minutes  from 
each  other)  and  performed  at  the  Acqua  Alta  Ocea¬ 
nographic  Tower  (AAOT)  in  the  northern  Adriatic 
Sea  during  clear  sky  conditions.  These  measure¬ 
ments  were  made  from  June  2002  to  November 
2008  in  a  variety  of  observation  conditions  indicated 
by  sun  zenith  0S  varying  in  the  range  of -23-70°  and 
chlorophyll-a  concentration  (as  quantified  through 
high  precision  liquid  chromatography)  in  the  range 
of  -0.2-3. 3  mg/m3. 

Above-water  data  were  produced  with  modified 
CE-318  (CIMEL,  France)  radiometers  [60]  integrated 
in  the  Ocean  Color  component  of  the  Aerosol  Ro- 
botic  Network  (AERONET-OC).  These  water-leaving 
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Fig.  6.  (Color  online)  (Top  panel)  Comparison  between  measured  and  derived  Lm(X,ft)/Lw(X,0°)  for  selected  wavelength  and  viewing 
angle  of  measurements  made  in  the  Mediterranean  Sea.  Red  triangle  represents  measurements  from  the  NuRADS;  blue  diamond  repre¬ 
sents  values  calculated  from  i?re  (IOPs-centered  approach).  (Bottom  panel)  Comparison  between  measured  and  derived  LW(X ,  Cl)/Lw(X ,  0°) 
for  n  in  the  remote  sensing  domain  (blue  circle,  IOPs-centered  approach,  green  square,  Chl-centered  approach). 
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radiance  data,  Lw(k,  6S  40°,  90°),  were  determined  at 
a  40°  viewing  angle  and  90°  relative  azimuth  with 
respect  to  the  sun  [61].  Water-leaving  radiance 
data  from  in-water  radiometric  measurements, 
LW(A,  0S  0°,  0°),  were  determined  from  subsurface  up- 
welling  radiance  Lu(k,0s  0°,0°)  extrapolated  from 
near-surface  nadir-view  measurements  performed 
with  the  Wire-Stabilized  Environmental  Profiling 
Radiometer  (WiSPER)  system  [62].  Spectral  differ¬ 
ences  between  the  two  datasets  were  minimized 
applying  a  band-shift  correction  relying  on  site- 
specific,  bio-optical  algorithms  [63]. 

R™(A,0s  40°,  90°)  was  calculated  for  each  set  of 
Lw(k, 0s ,40°,  90°)  using  theoretical  determinations 
of E<f(0+,  A, 0S)  for  clear  sky  conditions  [60],  and  then 
fed  to  the  QAA  (Section  4)  to  derive  a  and  b ^  values. 
With  these  a  and  b ^  values,  and  the  values  of  (G% (O), 
Gf(a)$  Gg(fl),  and  Gf(fl))  for  n  =  (0s,Oo,Oo)  and 
(1=  (0S,  40°,  90°),  0°,  0°)/i?rs(A,  40°,  90°) 

were  calculated  for  each  Lw(k ,  6S,  40°,  90°),  and  then 
Lw(A,0s,0*,0*)  were  determined  as  the  product  of 
R^Os.O^O^/R^  0S,4O°,  90°)  and  Lw(A,0s,  40°, 
90°). 

The  L^(A,^5,0°,0°)  calculated  from  the  AERO- 
NET-OC  Lw(k, 0^,40°,  90°)  are  compared  with  the 
WiSPER  Lw(A,0s,Oo,O°),  and  results  are  presented 
in  Fig.  8  (blue  symbols)  for  wavelengths  of  411, 
443,  490,  and  555  nm.  The  averages  of  signed  percent 
differences  (determined  to  quantify  biases  between 
compared  vales)  are  1.9%,  -0.8%,  -3.2%,  and 
-3.1%  at  the  four  wavelengths,  respectively.  The 
averages  of  the  absolute  of  percent  differences 


(determined  to  quantify  scattering  of  compared 
values)  are  7.8%,  6.9%,  8.0%,  and  6.8%.  If  the  BRDF 
correction  is  conducted  with  the  Chl-based  approach 
(green  symbols  in  Fig.  8,  where  Chi  was  estimated 
with  a  regional  algorithm  [64]),  the  averages  of 
signed  percent  differences  are  2.3%,  -0.0%,  -1.5%, 
and  -1.6%;  and  the  averages  of  the  absolute  of  the 
percent  differences  are  8.2%,  7.4%,  8.4%,  and  7.1%, 
respectively,  similar  as  that  from  the  IOP-based  ap¬ 
proach.  This  is  likely  due  to  the  relatively  small 
angular  effects  between  (40°,  90°)  and  (0°,  0°)  [1,5], 
as  without  a  BRDF  correction,  the  average  of  the 
signed  percent  differences  was  approximately  -2.2%, 
-5.3%,  -8  2%,  and  -8.2%,  respectively,  while  the 
absolute  of  percent  differences  are  ~8%  at  all 
wavelengths.  The  above  comparisons  indicate  that 
both  the  IOP-based  and  Chl-based  BRDF  corrections 
did  reduce  bias  but  not  appreciably  the  scattering 
among  compared  data  when  considering  Lw(k,0s, 
0°, 0°)  from  measured  LU,(A,^S,40°,90°).  The  latter 
insignificant  decrease  in  the  average  absolute 
percent  differences  is  explained,  at  least  in  part, 
by  uncertainties  characterizing  measurements  of 
water-leaving  radiance  through  two  fully  indepen¬ 
dent  instruments  and  methods. 

3.  Sensitivity  of  the  IOP-Based  Approach  on  the 
lOPs  Retrieval 

Because  of  various  assumptions,  approximations, 
and  empiricism,  a  full  set  of  IOPs  cannot  be  retrieved 
precisely  from  the  measurement  of  R^  (or  Lw).  To 
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Fig.  8.  (Color  online)  Comparison  between  measured  and  derived  L^^^O^O0)  for  measurements  made  at  the  AAOT.  The  x  axis 
represents  LW(A,0S<O° ,  0°)  measured  by  the  WiSPER,  while  they  axis  represents  Lw(jlf0s,OA,Oo)  calculated  from  Lw(X,&s,40°, 90°)  which 
was  measured  by  a  CE-318  radiometer  (blue  circle,  IOPs-centered  approach;  green  square,  Chl-centered  approach). 


understand  the  impact  of  lOP’s  uncertainties  on  the 
IOP-based  BRDF  correction,  for  example,  Fig.  9  pre¬ 
sents  the  change  of  Lw(Qs,Qv,(p)/Lw(6s,0° ,  0°)  corre¬ 
sponding  to  the  change  of  a(Ao)  for  the  two  sample 
measurements  in  Subsection  5.B.I.  We  show  this  be¬ 
cause  a(A0)  is  a  key  quantity  in  the  QAA  inversion 
scheme  for  the  sequential  calculation  of  other  optical 
properties.  Errors  or  uncertainties  of  the  other 
derived  properties  are  highly  dependent  on  the 
accuracy  of  a(A0)  but  not  uniformly  for  different 
waters  [47]. 

For  the  measurements  made  in  the  Mediterranean 
Sea  (blue  water),  a  ±10%  uncertainty  of  a(550)  was 
used  for  the  sensitivity  evaluation  (Fig.  9,  top  panel), 
while  a  ±20%  uncertainty  of  a  (550  )  was  used  for  the 
Monterey  Bay  water  (Fig.  9,  bottom  panel).  For  the 
two  waters  and  the  two  wavelengths  (412  and 
550  nm),  it  is  found  that  the  impact  of  10%  or  20% 
uncertainty  of  a(X0)  on  the  ratio  of  Lw(Qs,Qv,cp)/ 
Lw{Qs,  0°,0°)  is  limited.  Generally,  an  increase  of 
a(Ao)  will  increase  the  ratio  of  Lw (0§ ,Qv,(p)/ 
LW(0S,  0°,0°);  while  a  decrease  of  a(>lo)  will  decrease 
the  ratio  of  Lw(0s,0v,<p)/Lw{Qs,  0°,  0°)  in  similar 
magnitude.  For  the  blue  water,  a  10%  increase  of 
a(550)  resulted  in  2.4%,  2.7%,  3.2%,  and  3.8% 


increase  of  Lw(0s,Ov,<p)/Lw(0s,Oo,O°)  for  wave¬ 
lengths  of  412,  440,  490,  and  550  nm,  respectively, 
which  are  smaller  or  equivalent  to  the  measurement 
uncertainties  [65].  For  the  green  water,  although 
values  of  a (550)  were  altered  by  20%,  the  resulted 
impact  on  the  Lw(0Sy0Ln<p)/Lw(0sy  0°,0°)  ratio  is 
about  or  less  than  1%  for  412,  440,  490,  and  550  nm, 
respectively. 

The  limited  sensitivity  of  the  Lw(0g10v1<p)/ 
LW(6S,  0°,0°)  ratio  to  the  uncertainty  of  the  IOPs  is 
because  (1)  R re  is  generally  a  function  of  bb/(a  ±  h b) 
(e.g.,  Eqs.  (6)  and  (14)).  When  we  use  an  model  to 
analytically  derive  a  and  bbl  although  the  derived  a 
or  bb  values  may  be  inaccurate  in  some  degree  but 
the  retrieved  bb/{a  +  bb)  maintains  the  same  (e.g., 
the  QAA  scheme).  (2)  The  angular  variation  of 
RTS{0s,0v1(p)  is  primarily  governed  by  the  G(H)  coef¬ 
ficients  (Eq.  (14)),  which  is  determined  for  given 
particle  phase  functions. 

The  larger  uncertainty  in  the  Lw (6S )0v,(p)/ 
Lw(&Si  0°,0°)  ratio  of  the  blue  water  (comparing  to 
that  of  the  green  water)  is  due  to  a  lower  contribution 
of  particle  backscattering  to  the  upwelling  radiance 
with  respect  to  molecular  backscattering  (bbp/bb  was 
-25%  at  412  nm).  As  a  result,  the  relative  weighting 
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Fig.  9.  Scatter  plot  of  Lw(\,  Q)/Lw{\,  0°)  for  different  values  ofa(Ao).  The  x  axis  represents  the  Lw(k ,  Q)/Lw{\ ,  0°)  ratio  derived  with  a(550) 
in  the  current  QAA  scheme;  the  y-axis  represents  the  Lw(* l,  ft)/Lw(A,  0°)  ratio  derived  after  altering  the  estimated  a (550)  Top  panel,  the 
blue  water  in  the  Mediterranean  Sea,  a (550)  was  altered  by  ±10%;  bottom  panel,  the  green  water  in  the  Monterey  Bay,  a (550)  was  altered 
by  ±20% 


of  molecular  scattering  will  be  sensitive  to  the  de¬ 
rived  bbp  value.  Because  molecular  scattering  and 
particle  scattering  have  quite  different  angular 
shapes,  uncertainty  in  the  b^p  value  (resulted  from 
uncertainty  in  a(A0),  for  instance)  will  thus  affect 
the  ratio  of  Lw(0s,6lnq>) /LW(0S,  0°,0°),  which  is 
highly  dependent  on  the  shape  of  the  VSF.  For  the 
green  water,  however,  contribution  from  molecular 
backscattering  was  ~10%  of  the  total  backscattering, 
thus  the  phase  function  shape  of  molecular  scatter¬ 
ing  had  a  small  impact  on  the  shape  of  the  VSF, 
and  a  smaller  variation  in  the  Lw(6s,0v,q>)/ 
LW(0S,  0°,0°)  ratio  due  to  the  uncertainty  of  a(A0) 
was  found.  These  results  and  analyses  indicate 
slightly  larger  uncertainties  in  the  BRDF  correction 
in  oceanic  waters,  but  this  uncertainty  will  generally 
be  within  measurement  uncertainties  as  long  as  the 
particle  phase  function  is  well  defined. 

6.  Summary 

A  system  centered  on  IOPs  has  been  developed  for 
the  correction  of  angular  effects  in  water-leaving 


radiance  of  optically-deep  waters.  This  system  is 
not  restricted  to  Case  1  waters  and  focuses  on  the  de¬ 
rivation  of  the  total  absorption  and  backscattering 
coefficients  from  the  measured,  angular,  and  spectral 
remote-sensing  reflectance.  The  system  is  also  ap¬ 
plicable  to  multispectral  sensors  (e.g.,  SeaWiFS, 
MODIS,  and  MERIS)  and  hyperspectral  sensors. 
The  system  was  applied  to  both  a  numerically  simu¬ 
lated  dataset  and  a  field  measured  dataset  to  test 
and  evaluate  its  performance,  and  it  is  found  that  ap¬ 
plication  of  the  system  produced  adequate  estimates 
of  the  zenith  radiance  from  that  measured  at  other 
angles.  These  results  suggest  that  the  system  could 
be  applied  to  ocean  color  remote  sensing  to  (1)  derive 
IOPs  of  both  oceanic  and  coastal  waters  from  the 
angular  remote  sensing  reflectance  (or  angular 
water-leaving  radiance),  and  (2)  generate  normalized 
water-leaving  radiance  from  the  derived  IOPs  for 
calibration  and  validation  activities. 

HoweVer,  it  is  necessary  to  keep  in  mind  that  the 
accuracy  of  all  analytical  models  (including  the  mod¬ 
el  developed  and  used  here)  describing  the  remote 
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sensing  reflectance  is  limited  by  the  accuracy  of  the 
numerical  representation  of  the  radiative  transfer  in 
natural  waters.  More  importantly,  all  models  expli¬ 
citly  or  inexplicitly  include  an  assumption  about 
the  angular  shape  of  the  particle  scattering  function. 
This  angular  shape  plays  an  important  role  in 
forming  the  angular  distribution  of  the  upwelling 
radiance,  while  at  the  same  time  it  cannot  be  accu¬ 
rately  determined  yet  from  a  remote  sensing  mea¬ 
surement.  Because  the  exact  angular  shape  of 
particle  phase  function  is  not  known  a  priori  for  a 
given  water  body  in  remote  sensing,  and  the  model 
accuracy  is  no  better  than  the  accuracy  of  the  numer¬ 
ical  simulations,  correction  of  the  angular  variation 
of  water-leaving  radiance  based  on  the  IOPs  or  other 
properties  is,  to  the  best,  a  first  order  improvement. 
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